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Copper  metal  leaves  are  experimen tally  grown  in  two  dimensions  by  its 
electrodeposition  in  a  field  of  circular  symmetry  fran  a  solution  of  snail 
thickness.  The  determined  fractal  (D)  and  spectral  (d)  dimension  of  the 
deposited  aggregates  is  compared  with  that  determined  for  the  simulated  pattern 
calculated  ty  using  the  diffusion  limited  aggregation  method  of  Witten  and 
Sander  using  a  square  lattice.  Good  aggreement  is  found  for  the  fractal 
dimension  but  the  value  of  the  spectral  dimension,  "d,  is  found  to  be  40  %  larger 
for  the  observed  pattern.  This  results  from  the  fact  that  the  observed  pattern 
has  more  radial  symmetry  than  the  calculated  one.  This  implies  that  the 
chemical  length  the  center  to  the  perimeter  of  the  pattern  is  minimized 

during  the  growth.  Possible  causes  for  the  observed  deviation  are  discussed. 
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INTRODUCTION 

The  Witten  and  Sander  model  of  diffusion  limited  aggregates  (DLA)1  is  found 
to  be  useful  in  describing  the  dendritic  growth  phenomena.  The  aggregates  which 
are  grown  in  this  manner  are  fractal  ,  with  fractal  dimension  D  which  depends 
only  on  the  Euclidean  dimension  d  of  the  growth.  It  has  been  shown  that  this 
model  can  be  applied  under  certain  circumstances  not  only  to  the  solidification 
of  a  pure  or  a  binary  liquid^,  but  also  to  the  diffusion4  of  nonviscous  liquid 
into  a  viscous  one.  It  has  also  been  applied  to  describe  the  electrodeposition 
process'*.  Indeed,  Matsushita  et  al.®  have  demonstrated  that  electrodeposition 
of  zinc  metal  can  grew  aggregates,  the  fractal  dimension  of  which  corresponds  to 
the  value  expected  f ran  the  Witten  Sander  model. 

This  letter  is  aimed  at  reporting  experimental  results  on  copper 
electrodeposition  which  tends  to  confirm  the  generality  of  the  shape  of  the 
aggregates  grewn  by  this  method.  Hcwever,  it  turns  out  that  a  careful  analysis 
of  the  electrodeposited  pattern  demonstrates  that  this  aggregate  exhibits  a  more 
radial  symmetry  than  that  calculated  using  the  ELA  theoretical  model  on  the 
square  lattice.  We  interpret  this  feature  in  terms  of  a  process  which  tends  to 
minimize  the  chemical  length7  connecting  the  center  and  the  circumference  of  the 
pattern,  and  the  effect  of  the  geometric  lattice  model  (i.e. ,  square  lattice  for 
DLA)  where  the  aggregation  processes  are  carried  out. 

EXPERIMENTAL  SET  UP  AND  RESULTS 

A  brass  cathode  (5  mm)  diameter  with  an  emmersed  tip  of  2  irm  diameter 
located  at  the  center  of  a  circular  copper  anode  electrode  of  a  diameter  of 
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14  an  and  5  inn  high.  Both  electrodes  are  enmersed  in  an  electrolyte  CuSO^ 
solution.  The  whole  setup  is  contained  in  a  glass  container  having  a  flat  and 
horizontal  bottan.  The  tip  of  the  cathoda  confines  the  initiation  of  the  growth 
of  the  aggregate  f ran  the  exact  contact  point  between  the  glass  vessel  and  the 
tip  of  the  cathod.  The  centering  of  the  cathod  in  order  to  have  the  circular 
symmetry  was  made  by  hand  and  its  precision  was  better  than  1  mm.  The 
electrolyte  was  a  water  solution  of  cupric  sulfate  (CuS04) ,  the  concentration  of 
which  has  been  varied  from  0.1  to  1.0  M.  The  thickness  of  the  electrolyte  in 
the  vessel  has  been  also  varied  f  ran  3  mm  to  7  mm/  so  that  the  liquid  sanetimes 
covered  the  anode  and  sanetimes  did  not.  More  importantly/  we  expect  that  the 
aggregation  process  is  a  2  dimensional  process  for  length  scales  larger  than 
this  thickness  and  3  dimensional  one  on  the  contrary.  Neither  the  thickness  of 
the  electrolyte  nor  the  salt  concentration  is  found  to  have  observable  effect  cm 
the  shape  of  the  aggregate  at  large  length  scale.  The  electrical  potential 
difference  between  the  electrodes  was  kept  constant  during  any  growth  experiment 
and  was  varied  fran  2  V  to  15  V. 

At  the  beginning  of  the  experiments  we  have  checked  the  variation  of  the 
potential  inside  the  solution  and  measured  the  resistance  of  the  electrode  as  a 
function  of  the  salt  concentration/  of  the  applied  voltage  and  of  the  thickness 
of  the  solution.  These  dependences  were  normal.  Pictures  of  growing  aggregates 
have  been  recorded  during  the  growth  using  a  video  camera  placed  below  the 
vessel.  They  have  been  digitized  and  analyzed  using  a  computer  (VAX  11/780). 
The  analyzing  procedure  is  similar  to  that  used  in  the  zinc  case®. 
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The  picture  of  a  typical  copper  aggregate  is  given  in  Fig.  la  to  Id  for 
different  stages  of  the  growth.  The  observed  pattern  is  formed  f ran  branches  of 
characteristic  width,  RQ  of  the  order  of  1  -  5  rtm.  It  looks  very  much  like  the 
leaves  of  electrodeposited  zinc  observed  by  Matsushita  et  al.®  for  radii  larger 
than  Rq.  It  is  found  experimentally  that  Rq  depends  on  the  thickness  of  the 
electrolyte  and  on  the  applied  potential.  also  is  found  to  fluctuate  f ran 
pattern  to  pattern  so  that  no  predeterminated  pattern  can  be  obtained. 
Moreover,  slight  non-linearity  in  the  digitization  process  shows  relatively  more 
copper  in  the  oenter  than  it  is  observed.  Following  the  method  used  by  these 
authors  to  analyze  their  pattern®,  we  have  computed  the  variation  of  the 
integrated  dark  zone  in  the  pattern  as  a  function  of  the  distance  R  from  the 
electrode  center.  Assuming  that  the  integrated  intensity  of  the  dark  zone  at 
distance  R  is  proportional  to  the  mass  of  the  copper  at  that  distance,  a  best 
fit  procedure  lias  been  used  to  determine  the  fractal  dimension  of  the  aggregate: 

M(R)~rd  (1) 

in  the  range  of  interest  (R  >  RQ)  (cf .  Fig.  2) . 

The  average  value  of  D  obtained  for  the  different  aggregates  is  found  to  be 
1.6  0.2.  This  reproducibility  in  copper  deposition  is  poorer  than  the  5  % 
found  in  the  case  of  the  electrodeposition  of  Zinc®  and  is  in  agreanent  with 

O 

previous  results  .  It  might  result  from  the  fact  that  a  mixture  of  two 
substance  seems  to  be  deposited  in  the  case  of  copper,  most  likely  copper  and 
copper  oxide.  In  spite  of  this  difference  f ran  the  Zink  case,  D  seems  to  be 
similar  and  both  seem  to  shew  a  more  centrally  symmetric  pattern  than  that  was 
calculated  with  diffusion  limited  aggregation  model  (Fig.  le)  on  a  square 
lattice. 
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DISCUSSION 

We  will  assune  that:  1)  There  is  neither  a  charge  acamulation  in  the 
electrolyte  nor  a  potential  drop  at  the  interface;  and  2)  the  current  density 
at  the  interface  should  be  proportional  to  the  electric  field.  If  we  new  assume 
that  the  copper  leaves  are  perfect  conductors,  the  equation  governing  the  motion 
of  a  random  walker  in  the  Witten- Sander  model  can  be  used  to  describe  the  change 
in  the  electrical  potential  inside  the  electrolyte  during  the  growth  out  of  the 
electrolyte  solution.  This  model  yields  fractal  dimension  consistent  with  the 
experimentally  determined  value. 

However  our  experiments  have  detected  seme  difference  between  the  model  and 
the  experiment.  Before  giving  a  physical  meaning  to  this  difference,  it  is 
first  necessary  to  quantify  it.  This  requires  to  introduce  the  notions  of 
chemical  length13  Lc,  between  the  origin  and  a  point  at  a  distance  R.  This  is 
simply  the  length  of  the  actual  path  that  connects  this  point  to  the  origin. 
According  to  this  definition,  the  shorter  the  chemical  length  L  is,  the  less 
tortuous  the  path  is  and  more  radial  the  pattern  looks.  Measuring  the 
statistics  of  the  chemical  length  Lc(R)  as  function  of  the  Euclidean  distance  R 
between  the  center  of  the  pattern  to  another  point  for  both  the  experimental 
aggregate  and  the  simulated  one  diould  enable  us  to  quantify  their  difference  in 
the  following  manner.  It  is  shewn  that  L  (r)  obeys  a  power  law  R  ; 

L  (R)~R*  (2) 

v 

since  the  space  is  fractal,  o<  depends  on  the  geometry  of  the  fractal  and  can 
take  any  value  larger  than  1  (it  cannot  be  smaller  than  1  since  the  anallest 
distance  between  any  two  point  is  the  straight  line  which  lead  to  L  ~R) .  By 
measuring  o(  for  both  the  experimental  and  simulated  aggregates  we  may  quantify 


the  difference  between  than. 


In  order  to  relate  the  o(  exponent  to  other  exponent  already  defined  in  the 
case  of  a  fractal  space,  let  us  first  notice  that  in  the  computer  simulated  and 
the  experimental  patterns,  there  is  generally  only  one  path  which  connects  two 
points  (i.e. ,  there  is  no  circle  or  closed  loop  existed) ;  also  recall  that  in  a 

Q  =1  4 

fractal  space  the  spectral  dimension  ,  d,  governs  the  number  of  distinct  sites 
N(t)  that  a  random  walker  has  visited  fran  time  0  to  time  t: 

N(t)~t^2  (3) 

Under  these  circumstances  and  when  d  is  analler  than  2,  the  exploration  can  be 
considered  as  compact10  which  means  that  ary  point  separated  (i.e.,  frcrn  the 
center)  by  a  chemical  distance  L  snaller  than  1_  has  been  reached  by  the  random 
walker.  L  is  proportional  to  the  square  root  of  the  ellapsed  time  as  in  a 
classical  one  dimensional  diffusion  problem.  Combining  the  definitions  of  the 
fractal  dimension  (Eq.  1)  and  spectral  dimension  (Eq.  3)  together  with  this 
consideration  leads  to 

N(t)  -  t3''2  =  (LJ3~  R°  (4) 

Comparing  this  equation  to  the  definition  of  the  exponent  u  (Eq.  2)  leads  to: 
oi  =  D/d  (5) 

One  way  to  compute  the  exponent  o<  is  then  to  determine  d  through  a  random 
walk  problem  (Eq.  3) .  This  was  the  approach  we  have  used.  We  first  began  with 
the  ELA  model  in  order  to  check  the  accuracy  of  our  program.  We  get 
d  =  1.1±0.2  for  the  Witten  Sander  model.  This  result  is  in  agreement  with 
previous  determinations11.  The  large  fluctuation  we  get  are  mostly  due  to  the 
small  number  (1000)  of  sites  which  constitute  our  aggregate.  In  the  case  of  the 
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DLA  “  1,5 

o(  electrodeposition 

which  are  very  different  from  one  another.  This  result  oonfirms  that  the 
experimental  pattern  has  more  radial  symmetry  than  that  predicted  by  the  ELA 
model  calculated  on  a  square  lattice. 

There  are  several  reasons  that  might  account  for  the  difference  between  the 

calculated  and  the  observed  d.  The  ELA  model  is  an  ideal  model  for  describing 

electrodeposition.  It  does  not  consider  any  transport  properties  which  oould  be 

very  important  for  the  real  physical  systems.  For  example,  in  the 

electrodeposition  process,  the  Cu  aggregates  have  loose,  porous  structure 

instead  of  compact,  crystalline  structure.  This  might  increase  its  resistance. 

Since  the  resistance  within  the  aggregate  is  determined  ty  the  chemical  length 

Lc  of  the  aggregate  branches,  for  a  certain  radius  f ran  the  center  of  the 

patten-,  radial  branches  are  expected  to  have  lower  resistance  than  branched 

ones.  This  leads  to  higher  current  density,  and  thus  faster  rate  of  deposition 

at  the  end  of  radial  branches  than  at  the  end  of  more  tortuous  branches  having 

the  same  distance  f ran  the  center.  Other  possible  causes  for  the  more  radial 

appearance  of  the  observed  pattern  could  simply  be  an  inhomogeneous  depletion  of 
+2 

the  Cu  concentration  during  the  deposition  process.  Due  to  the  fragile  nature 
of  the  two  dimensional  pattern,  stirring  was  not  carried  out. 

While  we  have  hlamed  the  disagreement  between  the  observed  and  calculated 
spectral  dimension  on  deviation  of  the  deposition  process  f  ran  the  ELA  model, 
sane  of  the  deviation  might  actually  result  f  ran  the  use  of  a  square  lattice  in 
the  simulation  calculation  of  the  pattern  according  to  the  ELA  model.  It  has 
already  been  pointed  out^  that  the  value  of  the  fractal  dimension,  D,  depends 
on  the  type  of  lattice  used  in  these  types  of  calculations.  It  might  be  that 
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the  value  of  the  spectral  dimension  is  also  sensitive  to  the  type  of  the  lattice 
used. 

Studies  are  new  in  progress  and  aimed  at  examining  the  sensitivity  of  the 
observed  fractal  and  spectral  dimension  to  variation  in  the  different 
experimental  conditions.  We  are  also  examining  the  sensitivity  of  the  spectral 
dimension  of  the  calculated  pattern  on  the  type  of  lattice  used  in  the 
calculation. 
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Figure  Caption 


Fig.  Is  Digitized  form  of  the  observed  elect  rode  posited  copper  aggregate 
at  four  different  stages  of  the  growth:  voltage  used  is  2.0  volt,  copper  ion 
concentration  is  0.5  M,  the  complete  time  of  growth  is  50  minutes.  Fig.  le  is 
the  simulation  of  a  aggregate  grown  by  diffusion  limited  aggregation  model.  The 
first  four  patterns  exhibit  a  more  radial  symmetry  than  the  simulated  one. 


Fig.  2:  A  Plot  of  In  (N(R) )  vs.  In  R  for  different  observed 
electrodeposited  copper  patterns.  It  is  assumed  that  N(R)  is  proportional  to 
the  area  of  the  dark  area  of  the  digitized  form  of  the  pattern  in  a  circle  of 
radius  R.  From  the  slope  of  the  line,  the  fractal  dimension  D  is  determined. 


Fig.  3  s  The  determination  of  the  spectral  dimension  for  the  observed  (a) 
and  the  calculated  (b)  patterns. 
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